Indication of the ferromagnetic instability in a dilute two-dimensional electron system 
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The magnetic field Be, in which the electrons become fully spin-polarized, is found to be pro- 
portional to the deviation of the electron density from the zero-field metal-insulator transition in a 
two-dimensional electron system in silicon. The tendency of B^ to vanish at a finite electron density 
suggests a ferromagnetic instability in this strongly correlated electron system. 
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At sufficiently low electron densities, an ideal two- 
dimensional (2D) electron system becomes strongly cor- 
related, because the kinetic energy is overpowered by en- 
ergy of electron-electron interactions (exchange and cor- 
relation energy). The interaction strength is normally 
described by the Wigner-Seitz radius, r^, which is equal 
to the ratio of the Coulomb and the Fermi energies as cal- 
culated for 2D band electrons. In the low-density limit, 
at Ts ^ 35, the ground state of the system is expected 
to be a Wigner crystal characterized by spatial and spin 
ordering 0|. The spin ordering may survive at higher 
electron densities (lower r^) up to the threshold deter- 
mined by ferromagnetic (Stoner) instability In this 
strongly- interacting limit {rg >> 1), all results given by 
theoretical approaches are very approximate, not to men- 
tion that in real 2D electron systems, the influence of dis- 
order needs to be taken into account, which complicates 
the problem drastically. In strongly disordered 2D sys- 
tems, the influence of the disorder can dominate the in- 
teraction effects leading to a disorder-driven localization, 
whereas in the least disordered 2D systems, the metal- 
insulator transition may be caused by interaction effects 

In magnetic fields parallel to the 2D electron plane, the 
spin effects should dominate the properties of a 2D elec- 
tron system once the orbital quantization is quenched. 
Indeed, the resistance of a 2D electron gas in silicon 
metal-oxide-semiconductor field-effect transistors (MOS- 
FETs) was found to be isotropic with respect to in-plane 
magnetic field, B, and rise steeply with the field saturat- 
ing to a constant value above a critical magnetic field 
Be which depends on electron density jj]. Moreover, 
an analysis of Shubnikov-de Haas oscillations in tilted 
magnetic fields has established recently that the field Be 
corresponds to the onset of full spin polarization of the 
electron system 

In this paper, we study low-temperature parallel-field 
magnetotransport in a 2D electron system in silicon in a 



wide range of electron densities, Ug. We find that the sat- 
uration (or polarization) magnetic field. Be, is a strictly 
linear function of Us'- Be (x (jig — ric) where Ue is the 
critical electron density for the B — metal-insulator 
transition. Since above ne, the 2D bandtail of localized 
electrons is negligibly small, as inferred from the low- 
field/low-temperature Hall effect measurements, vanish- 
ing Be at a finite electron density gives evidence in favor 
of the existence of a ferromagnetic transition (at ris — rie) 
in this strongly correlated 2D electron system. 

Measurements were made in a rotator-equipped Ox- 
ford dilution refrigerator with a base temperature of 
K, 30 mK on high-mobility (lOO)-silicon samples similar 
to those previously used in Ref. 0. The resistance was 
measured by a standard 4-terminal low-frequency tech- 
nique. Excitation current was kept low enough to ensure 
that measurements were taken in the linear regime of 
response. Contact resistances in our samples were mini- 
mized by using a split-gate technique that allows one to 
maintain a high electron density in the vicinity of the 
contacts regardless of its value in the main part of the 
sample. In this paper we show results obtained on a 
sample with a peak mobility close to 3 m^/Vs at 0.1 K. 

Typical curves of the low-temperature magnetoresis- 
tance p{B) in a parallel magnetic field are displayed in 
Fig. ^ The resistivity increases approximately quadrat- 
ically with field until it saturates at a constant value 
above a certain density-dependent magnetic field. In 
the vicinity of the metal-insulator transition, the magne- 
toresistance strongly depends on temperature T, as was 
reported, e.g., in Ref. As one moves away from the 
transition, however, this temperature dependence weak- 
ens and eventually disappears at very low temperatures 
where the metallic resistivity at B = saturates and 
becomes independent of temperature. Therefore, as Us 
is lowered, the low-temperature limit is realized for pro- 
gressively narrower initial intervals on the curve p{B). 
The data discussed in this paper are obtained in this 
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FIG. 1. Low-temperature magnetoresistance in parallel magnetic 
fields at difli'erent electron densities above the critical density for 
the B = metal-insulator transition. The lowest density curve is 
outside the scaling region as described in the text. 



FIG. 2. Scaled curves of the normalized magnetoresistance at dif- 
ferent Ug vs B/Bq. The electron densities are indicated in units 
of 10^^ ia~^ . Also shown by a dashed line is the normalized mag- 
netoresistance calculated in Rcf. [hd. 



low-temperature limit where the magnetoresistance be- 
comes temperature-independent. 

In Fig. H, we show how the normalized magnetoresis- 
tance, measured at different electron densities, collapses 
onto a single curve when plotted as a function of B / Be- 
The scaling parameter. Be, has been normalized to cor- 
respond to the magnetic field at which the magnetore- 
sistance saturates (within the accuracy with which the 
latter can be determined) and hence full spin polariza- 
tion of the electrons is reached [^J^ . The observed scal- 
ing is remarkably good for B/Bc < 0.7 in the electron 
density range 1.08 x 10^^ — 10^^ m^^, although with in- 
creasing Us the scaled experimental data occupy progres- 
sively shorter intervals on the resulting curve. Both at 
B/Bc > 0.7 and outside the indicated range of elec- 
tron densities, the scaled data start to noticeably de- 
viate from the universal curve. In particular, the scal- 
ing breaks down when one approaches {rig < l.Sric) the 
metal-insulator transition which in this sample occurs 
at zero magnetic field at ric = 8 x lO^'' m~^ [^. This 
is not surprising as the magnetoresistance near ric de- 
pends strongly on temperature, see above. We note that 
the observed scaling dependence is described reasonably 
well by the theoretical dependence of p/p{0) on the de- 
gree of spin polarization ^ = gmfiBB/nh^ns = B/Bc 
(where both the g factor and the effective mass m may 
be enhanced due to interactions, see below) as predicted 
by the recent theory |l2| based on the spin-polarization- 
dependent screening of a random potential deeply in the 
metallic regime, see Fig. ||. 

In Fig. 1^, we show the first important and unexpected 
result of this paper: with high precision. Be is propor- 
tional to the deviation of the electron density from its 
critical value, i.e., to (n^ — ric), over a wide range of 
electron densities. In other words, the field, at which the 
magnetoresistance saturates, tends to vanish at ric- We 
emphasize that our procedure provides high accuracy for 
determining the behavior of the field of saturation with 



electron density, i.e., the functional form of Bc{ns), even 
though the absolute value of Be is determined not so 
accurately. Note that at Ug above 2.4 x 10^^ m^^, the 
saturation of the resistance is not reached in our mag- 
netic field range; still, the high precision of the collapse 
of the high-density experimental curves onto the same 
scaling curve as the low-density data allows us to draw 
conclusion about the validity of the obtained law Bcijig) 
over a much wider range of electron densities. 

The fact that the saturation (or polarization) field 
tends to vanish at Uc would be trivial if the density of 
the delocalized electrons approached zero at the metal- 
insulator transition. This crucial point can be verified 
with the help of the Hall effect measurements in the limit 
of weak magnetic fields and low temperatures |T^. As 
seen in the inset to Fig. ^, above Uc, the electron density 
obtained from the Hall effect measurements agrees well 
with the one determined from Shubnikov-de Haas oscil- 
lations, which shows that in the metallic phase, the 2D 
bandtail of localized electron states is small, i.e., all the 
electrons are delocalized [0. Then, the observed ten- 
dency for Be to vanish at a finite electron density gives 
a sign of the long-awaited ferromagnetic instability in 
strongly correlated 2D electron systems jTsj. 

Subject to the occurrence of a ferromagnetic (Stoner) 
instability at = Uc, it is easy to calculate the be- 
havior of Be vs Hs based on the functional form of the 
exchange and correlation energy which is valid in the 
strong-interaction limit. Apparently, the Fermi energies 
of a spin polarized and unpolarized 2D electron system 
with both interaction and Zeeman energy taken into ac- 
count must be equal. Writing the exchange and corre- 
lation energy in the explicit form, we get for the energy 
difference at the onset of full spin polarization and at 
B = 

^-7--/ ^-90PBBe = 0, (1) 
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FIG. 3. Dependence of the field _Bc on electron density. The dashed 
line is a linear fit which extrapolates to the critical electron density 
for the B = metal-insulator transition. The fit using Eq. (1) 
with ^ = 0.095 is shown by a solid line. A close-up view of the 
region near ric is displayed in the inset. 

where h is the Planck constant, the band mass nib — 
0.19 mo (mo is the free electron mass), e is the electron 
charge, k is the dielectric constant, the Lande g factor in 
bulk silicon is equal to go = 2, and /is is the Bohr mag- 
neton. The first term in Eq. (1) is given by the difference 
of the bare Fermi energies which differ by a factor of two 
originating from the spin degeneracy (valley degeneracy 
being taken into account). The interaction (second) term 
contains the unknown numerical factor 7 which is posi- 
tive since the exchange and correlation energy is always 
negative and its absolute value enhances with spin po- 
larization. Note that the interactions are expected to 
lead to corrections to both nih and go; the problem of 
the interaction-enhanced effective mass m and g factor 
lacks a definite answer so far. Note also that the dras- 
tic enhancement of the product gm (Fig. ^ determined 
from the data in Fig. ||, based on the above form of the 
spin polarization parameter ^ is consistent with the gm 
value obtained from the analysis of Shubnikov-de Haas 
oscillations in dilute sihcon MOSFETs 

As seen from Fig. |^, the resulting dependence Bdris) 
with the factor 7 = 0.095 determined from the condi- 
tion Bc{nc) = describes the experimental finding well 
enough. Formally, the best description of the experimen- 
tal data might be obtained with a bit smaller 7 = 0.090, 
which corresponds to the expected ferromagnetic transi- 
tion point just below ric, at = 7 x 10^^ m~^. Although 
the accuracy of the Be normalization hardly allows one 
to do so, the conclusions of the paper would still hold for 
that case. We emphasize, however, that the Bc{ns) de- 
pendence given by Eq. (1) is weakly superlinear, whereas 
our experiments yield strictly linear dependence. Clearly, 
theoretical efforts are needed to explore a possible rela- 
tionship between the Stoner instability and the metal- 
insulator transition in 2D. 

In summary, we have studied the low-temperature 
magnetoresistance in parallel magnetic fields in a wide 



FIG. 4. The product gm/gQlTlh vs Ug obtained from the data in 
Fig. ^, based on the above form of the spin polarization parameter 
^. The dashed line is a guide to the eye. The inset compares the 
electron densities determined from Shubnikov-de Haas oscillation 
and weak-field low-temperature Hall measurements. For the latter, 
averages were taken for opposite magnetic field directions. 

range of electron densities above the critical electron den- 
sity for the zero-field metal-insulator transition, ric- The 
normalized magnetoresistance is found to scale with Us 
defining the scaling parameter, Be, that corresponds to 
the magnetic field in which the full spin polarization is 
achieved. Over a wide range of electron densities, this 
scaling parameter changes in precise accordance with 
the relation B^ oc (ris — ric). Although in the metallic 
regime, the bandtail of the localized electrons has been 
found negligibly small, B^. tends to vanish at a finite elec- 
tron density. This gives evidence in favor of the existence 
of a ferromagnetic transition in this strongly correlated 
2D electron system. 

Recently, a conclusion about possible ferromagnetism 
in a dilute 2D electron system has been made by another 
experimental group [|T]. 
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